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ABSTRACT 



Dynamical studies of MHD turbulence on the one hand, and arguments 
based upon magnetic helicity on the other, have yielded seemingly contradictory 
estimates for the a parameter in turbulent dynamo theory. Here we show, 
with direct numerical simulation of three-dimensional magnetohydrodynamic 
turbulence with a mean magnetic field, B, that the constraint on the dynamo 
a-effect set by the magnetic helicity is time-dependent. A time-scale t c is 
introduced such that for t < t c , the a-coefficient calculated from the simulation is 
close to the result of Pouquet et al. and Field etal. , - ^ ( ( v • V x v) - (b • V x b) ) ; 
for t > t C) the classical result of the a-coefficient given by the Mean-Field 
Electrodynamics is reduced by a factor of 1/ (i? m |B| 2 /t>;! ms ) , as argued by 
Gruzinov & Diamond, Seehafer and Cattaneo & Hughes. Here, R m is the 
magnetic Reynolds number, v rms the rms velocity of the turbulence, r cor the 
correlation time of the turbulence, and B is in velocity unit. The applicability 
of and connection between different models of dynamo theory are also discussed. 
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1. Introduction 

The generation and amplification of magnetic field in many astrophysical systems 
are often attributed to the turbulent dynamo effect. Since the seminal papers by Parker 
(1955) and Steenbeck, Krause & Radler (1966), a whole body of theory, namely, Mean 
Field Electrodynamics (MFE) has been developed to explain the dynamics of magnetic field 
generation by helical turbulence in a conducting fluid (Moffatt 1978, or Krause & Radler 
1980). In such a fluid, the velocity field, v, stretches the magnetic field, b, in such a way 
that the correlation of v and b results in an electromotive force, (v x b) that amplifies the 
mean (large-scale) magnetic field, B, through the relation 

(v x b) = aB. (1) 

Here, the coefficient a represents the so-called a-effect (Moffatt 1978); it is calculated in 
MFE as olmfe ~ — ■ V x v), where r cor is the correlation time of the turbulence. 
MFE is a kinematic theory, in that the velocity field is prescribed and no back reaction 
of the magnetic field on the velocity field is considered. Therefore, its applicability 
to circumstances where the velocity field is affected by the growing magnetic field is 
questionable. Several authors have extended MFE to include the quenching of the a-effect 
due to the back reaction of the magnetic field. By numerically solving the spectral MHD 
equations using a closure method known as the EDQNM (Eddy-Damped Quasi-Normal 
Markovian) approximation, Pouquet, Frisch & Leorat (1976, hereafter PFL) find that 
a = j k akdk, where is determined by, in Fourier space, the difference between the spectra 
of the kinetic helicity correlation function, ( v ■ V x v) , and the current helicity correlation 
function, (b ■ V x b). A similar result was found by Field, Blackman & Chou (1999, 
hereafter FBC), who consider the back reaction of magnetic field by treating v and b on an 
equal footing, and give a model in which the a-coefficient (for later purposes, we call it a±) 
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can be expressed as 



«i = -j((v-Vxv)-(b-Vxb)) (2) 



for relatively small values of [B/v rm ^j 



2 



However, the nonlinear nature of the problem introduces so much difficulty that the 
effect of the back reaction is still under debate. One of the objections to the application of 
cci to the galactic dynamo has its root in the problem of large magnetic Reynolds number 
and small-scale fields. Vainshtein & Cattaneo (1992, see also Cattaneo & Vainshtein 1991) 
argue that for systems with large R m , a is reduced by a factor of R m from its kinematic 
value, i.e., 

a oc -j-, (3) 

as the small-scale magnetic field grows quickly to turn off the generation of magnetic flux 
at large scales. These authors applied the conservation of the square vector potential of 
2D ideal MHD. The conservation of magnetic helicity of 3D MHD in steady state was later 
studied by Seehafer (1994, 1995) who provided another model of the dynamo ct-effect in 
which the a-coefficient (for later purposes, we call it a 2 ) is given by 

a 2 = '—= 2 '-. (4) 

t> 

Because R m = Lv rms /\, a 2 oc 1/R m also. Because R m ^> 1 in astrophysical systems such 
as the Galaxy, both relation (3) and relation (4) suggest strong suppression of the dynamo 
a-effect. 

Another model of the dynamo a-effect was introduced by Gruzinov & Diamond (1994, 
1995, 1996) by studying how the conservation of magnetic helicity affects the result like 
relation (2). They realized that because Ohmic dissipation of the current helicity (b-Vxb) 
in (2) changes the magnetic helicity, the dynamics of the latter will affect a 1: which should 
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be modified to (for later purposes, we call it a 3 ) 

_ OLMFE /_x 

a3 ~ l + R m (B/v rms y 1 j 

Since R m is usually very large in astrophysics, this again implies that a, in steady state, is 
far smaller than its classical value, so that it is too small to be important for the generation 
of large scale magnetic fields. Note that in the limit of R m ^> 1 and v rms ~ b rms , a 3 is 
reduced to a 2 . 

a 2 and a 3 are different from ai in the following two aspects: first, a.\ does not 
explicitly show strong suppression of the a-effect for large R m ; second, neither of the 
dynamical studies of PFL and FBC, which give aii, explicitly considered the conservation 
of magnetic helicity that led to a 2 and a 3 . Under what conditions the magnetic helicity 
constraint, which is essential in the models of a 2 and a 3 , enters the derivation of a\ is still 
an open question. In other such words, if such magnetic helicity constraint can be relaxed 
in systems of large R m , which of these three models remains valid? To provide insight 
into this problem, in the following we study the dynamics of magnetic helicity using the 
magnetic helicity conservation equations, and the time dependence of a using a numerical 
simulation of MHD turbulence. We find that both the magnetic helicity development and 
the a-coefficient are time dependent. We find that the classical MFE result is valid up 
to a critical time that we calculate, and the magnetic helicity constrained result is valid 
thereafter. Which value to apply therefore depends on the circumstances. 

This paper is presented in the following structure: in Section 2 we provide a model 
for the time-dependence of the magnetic helicity development, and introduce a critical 
time, t c , to separate the two important stages of the magnetic helicity evolution; in Section 
3.1 we present our numerical model that is used to study the time-dependence of the 
dynamo a-effect; the numerical results are given in Section 3.2; in Section 4 we discuss 
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the implications of our analytic model and the numerical results, which are applied to the 
Galactic dynamo; conclusions are given in Section 5. 



2. The Time Dependence of Magnetic Helicity Dynamics 

We start from studying the dynamics of the magnetic helicity and its constraint on the 
dynamo a-effect. For a 3D incompressible MHD system, we separate the vector potential 
A into a large-scale part A and a small-scale part a. Similarly, we write the magnetic field 
as B = B + b and the velocity field as V = V + v. By un-curling the induction equation 
for B, 

r)B 

°— = V x (V x B) + AV 2 B, (6) 



we have the equation for A, 



dA 

— = V x B + AV 2 A - W- (7) 

at 



Here A is the magnetic diffusivity and ip is the scalar potential. Dotting (6) with a and 
(7) with b and summing the resulting equations together, we have the equation for the 
ensemble average of the small-scale magnetic helicity 

J^(a-b) = -2A(b-Vxb)-2(vxb)-B 

+ (V- (-b^'-ax (vxb + Vxb + vxB))), (8) 

where ip' is the fluctuating component of the scalar potential. The third term in (8) comes 
from the v x B term in both (6) and (7), which represents the interaction between the 
small-scale velocity field v and the large-scale magnetic field B. The physics of this term 
can be explained as follows. The line stretching, twisting and folding of B will produce 
b and a, therefore affect the generation and diffusion of (a ■ b) within a volume V. The 
divergence term shows that flux of magnetic helicity of certain sign can escape from the 



system through open boundaries (Blackman & Field, 2000). The equation for the large-scale 
magnetic helicity, A • B, can be also derived. The equations for B and A are 

<9B 



V x (V x B) + V x (v x b) + AV 2 B, (9) 
Vx B+ (vx b) + AV 2 A- V^. (10) 



dt 

dA „ _ 



dt 

Following similar procedures that led us to get (8), we have 

d _ _ 

A B = — 2AB ■ V x B + 2(v x b) • B 

C6 

+ (V- (-5$- A x (V x B+ (v x b)))). (11) 

Note that the third term in (11) has the opposite sign of the third term in (8), showing that 
the electromotive force ( v x b) generates A • B of the opposite sign of (a • b) . If we define 
the following two flux terms 

f = (b^' + ax (v x b + Vx b + v x B)), (12) 

F = (B0 + A x (V x B + (v x b))), (13) 
with (12) and (13), we may re-write (8) and (11) in the limit of A — > in the form 

A (a • b) = ^(a • b) + V • f = -25 • B, (14) 

d _ 

D t (A • B) = —A • B + V • F = 2£ • B, (15) 
dt 

where S = (v x b) is the electromotive force. The conservation of total magnetic helicity, 
H m = (a • b) + A ■ B, immediately follows from the above two equations and reads 

V t H m = A(a-b)+A(A-B) = 0, (16) 

which holds for ideal MHD. 
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Seehafer (1994, 1995) related the dynamics of small-scale magnetic helicity, i.e., 
equation (8), to the quenching of the dynamo a-effect. By assuming stationarity of the 
MHD turbulence in a closed system, he argued that d t (a. • b) =0 and V • f = V • F = 0; 
therefore, by neglecting the time-dependent term and the boundary term in equation (8) 
and applying relation (1), he obtained relation (4) for a 2 . The same assumptions about 
the stationarity and the closedness of the system were made by Gruzinov & Diamond 
(1994, 1995, 1996), where they provide a modified a-coeflicient in the form of a 3 . a 3 can 
be regarded as an interpolation between the MFE result cimfe and the quenching result 
a = rni FE \2 ( see a l so Cattaneo & Vainshtein, 1991 and Vainshtein & Cattaneo, 1992), 
thus it is supposed to be valid for both large and small R m . The numerical simulation by 
Cattaneo & Hughes (1996) with a particular R m = 100 and different values of B supports 
this result. 

«2 and «3 are different from a± in that both of these two models show strong suppression 
of a for large R m . The large magnetic Reynolds number R m in real astrophysical systems 
constrains the dynamo a-effect to such a degree that, according to these two models, a 
would be too small to be important. However, from the above derivation of both a 2 and 0:3, 
one can see that three conditions must be met in order for such constraint to be effective in 
real astrophysical systems: 

• a large magnetic Reynolds number, 

• the system is in stationary state, and 

• the system is closed, i.e, no net flux of magnetic helicity flowing through the system. 

Because the numerical simulation of Cattaneo & Hughes (1996) satisfies all these conditions, 
it is interesting that their numerical results confirm a 3 in relation (5). 
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For real astrophysical systems, these three conditions may not be satisfied 
simultaneously. For example, Blackman & Field (2000) have argued that most astrophysical 
objects are open systems, and magnetic helicity can flow through the boundaries. For 
systems where B cannot be assumed constant, Bhattacharjee & Yuan (1995) suggested that 
a — > =^V ■ (k 2 V=^) for R m — > oo, where J = V x B and k 2 is a positive functional of the 
statistical properties of the MHD turbulence. 

Another assumption made in deriving a 2 and 0:3, the stationarity of the MHD 
turbulence, may also not be valid when the a-coeflicient in relations (4) or (5) is 
applied to real astrophysical systems. One scenario that we can imagine is that when 
impulsive, transient phenomena such as solar flares or supernova explosions happen, 
steady astrophysical systems will be disturbed. This means that, to correctly understand 
the relation between the dynamics of magnetic helicity and the dynamo ct-effect and its 
application to real astrophysical systems, one should study not only the stationary state 
where the magnetic helicity has been built up in the system, but also the non-stationary 
state when there is net magnetic helicity being built up in or thrown away from the system. 
During such non-stationary state of the magnetic helicity development, the velocity field has 
the freedom to stay independent from the growing magnetic field, and if so, the derivation 
of a 2 from just the induction equation may not give a complete picture for the dynamo 
a-effect in real astrophysical systems (Kulsrud 1999). On the other hand, both PFL and 
FBC, who give a± in dynamical studies, incorporate not only the momentum equation but 
also the induction equation into their models. Because a 3 can be derived by relating the 
(b • V x b) terms in relations (4) and (2), it is expected that 0:3 should include not only 
the dynamical considerations of PFL and FBC but also the magnetic helicity constraint of 
Seehafer. However, in order to derive a 3 in this way, one has to assume that cei is equivalent 
to «2 in order to get 0:3 in the form of relation (5). This means that the assumptions of 
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closedness, stationarity and large R m must be satisfied in the model of a 3 . 

If the stationarity assumption is removed, the magnetic helicity constraint shown in 
the models of a 2 and a 3 can then be relaxed, even though in these models the magnetic 
Reynolds number is large and the system is closed. And if the magnetic helicity constraint 
is relaxed, the model of a\ may give a complete picture of the dynamo a-effect. If the 
relaxation of the magnetic helicity constraint lasts for a limited period of time, the model 
of aii must be extended so that the magnetic helicity conservation is taken into account by 
introducing a 2 and a 3 . Note that a 2 ~ 0:3 for large R m , and we are only interested in large 
R m (or small A, as in astrophysics) case. 

To test our hypothesis of the non-stationary state of a 3D MHD system, we consider 
a closed system and neglect the boundary effects due to V • f and V • F. The resulting 
equations for the small and large scale magnetic helicity can be written as 



where relation (1) is used. The interpretation of equations (17) and (18) is that: the 
production of positive/negative large-scale magnetic helicity, A • B, by the dynamo a-effect 
is due to the production of negative/positive small-scale magnetic helicity, (a - b) and to 
dissipation, A(B- V x B). For systems with A — > 0, the total magnetic helicity, A-B+ (a-b), 
is conserved, and dynamo a-effect can be understood as a pumping effect that transfers 
small-scale magnetic helicity to large scales without generating any total magnetic helicity 
in the (closed) system. To illustrate such a pumping effect, consider the case that at certain 
time the small-scale magnetic helicity is zero, (a ■ b) = 0. A dynamo a-effect with positive 
a coefficient will pump positive magnetic helicity from small scales to large scales, and leave 




d 



(17) 



and 




(18) 
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negative magnetic helicity at small scales, so that after some critical time t c , (a • b) < 0, 
A • B > 0, and (a • b) + A ■ B = 0. 

To quantify these effects, consider the case in (17) that A is very small so that the 
diffusion term can be neglected. If at t — 0, (a - b) = 0, the integral of the resulting equation 
givesQ 

(a • b) (t) = -2B 2 I a(r)dr. (19) 
Jo 

For small times, when the constraint of magnetic helicity on the dynamo effect is not yet 
effective, we assume that 

a{t) = an, (20) 
the dynamical value calculated by the model of ai, so 

(a- b)(t) = -2B 2 a x t. (21) 

This should be valid for small values of t. However, as t increases, (a • b) (t) approaches 
the maximum that can be associated with the given small-scale magnetic energy It can be 
shown from the realizability condition (Moffatt 1978) that 

|(a-b)(t)| <2Lb 2 rms , (22) 

where L is the outer scale of the turbulence, and b rms is the rms small-scale field strength. 
In light of this, (20) is valid only up to a critical time t c given by 

t c = (23) 
ai B 



^ote that here we are assuming that B is a independent of time, as in our numerical 
simulation. 
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Since a x ~ 0(v rms ), where f rms is the rms velocity, and L/v rms — t eddy , the eddy turnover 
time, 

t c ^t eddy {^ . (24) 

A more realistic estimate of t c should include a correction factor, which modifies the above 
formula to t c = Ct eddy (^§p) . The correction factor C depends on the ratio v rms /a\ and 
the ratio |(a • b)\/2Lb 2 . ms . C can be estimated from numerical simulation results, as we will 
do in the next section. 

When t > t c , because there is no further small-scale magnetic helicity to draw on, 
d t (a ■ b) — > 0, and from (17) we have 

oB 2 = -A(b- V x b), (25) 

which implies that a no longer equals cki, but rather is determined by the fact that the 
only available source for the large-scale magnetic helicity being pumped by a is that which 
is being dissipated at small scales with appropriate sign. That is, a equals «2, hence is 
suppressed by the large R m (or small A) for the stage t > t c . To estimate the value of 
the suppressed a, a sp , we replace the current helicity of the small-scale magnetic field, 
(b • V x b), with the value corresponding to a maximally helical small-scale field (which we 
assumed to result from (20) at t — t c ), which is 

|(b.Vxb)|-L-\ 2 ms . (26) 

Hence 

(27) 

Since |ai| ~ v rms , and R m = Lv rms /\, this can be written as 

(28) 



1 (h \ 2 

1 / Orm.fi \ 



Rm \ B 
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This formula is consistent with a 3 when R m is large, and is identical to a 2 , the model 
by Seehafer (1994, 1995). However, notice that unlike previous derivations, ours depends 
critically on time; so 



Note that neither a 2 nor a 3 appears in relation (29) because they should be valid only 
for t > t c , when the magnetic helicity constraint takes effect. 

In the above discussion, we made the following assumptions that require further 
justification: (1) the system has an amount of small-scale magnetic helicity that can be 
estimated as Lbf ms , where L is the scale where the magnetic helicity is concentrated; (2) 
current helicity can be estimated as b^. ms /L; (3) v^ ms ~ ^ ms . To test these assumptions and 
calculate t c as a function of t eddy b1 ms / B , we apply direct numerical simulation to a 3D 
incompressible MHD system with periodic boundary conditions. The numerical model and 
simulation results are presented in next section. 

3. Numerical Model and Simulation Results 

3.1. Numerical Model 

Under an external force J 7 , the undimensionalized incompressible MHD equations can 
be written as (with Einstein summation convention) 




(29) 




(30) 




(31) 



diVi = djbj = 0, 



(32) 
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where v and A are the molecular viscosity and magnetic diffusivity, respectively. Note that 
we have written b in units of y/Aixp after dividing both sides of the Navier-Stokes equation 
and the induction equation by density p. If we use a hat, A, to denote discrete Fourier 
transform, and (g) to denote convolution, the above equations in Fourier space are 

(d t + uk 2 ^ vj = Pji \k m (-fy ® v m + b t ® S m ) + f x , (33) 

(d t + \k 2 ) bj = Pjt [ik m (vi <g> b m - v m <g> 6,)] , (34) 
km^m = kjbj = 0. (35) 

Here P is the projection operator defined as Pji = 5ji — In our simulation, we treat the 
system as a cube [0, 2n) x [0, 2n) x [0, 2n). The Cartesian coordinate of a grid point can be 
written as xi = ^1, y m = j^m, z n = ^n, for I, m, n — 0, 1, 2, N — 1. A point in Fourier 
space has coordinates k s = s,k p = p, k q = q, for s,p, q = — y, — ^f^, y — 1- Because we 
assume periodic boundary conditions, the surface terms V • / and V • F in equations (14) 
and (15) both vanish. 

Equations (33), (34) and (35) are numerically solved with the standard Fourier spectral 
method. Equations (33) and (34) are treated as ordinary differential equations for v and 
b. With the projection operator P, the divergence free condition (35) will be satisfied for 
t > as long as v and b are divergence free at t — 0. All our simulation runs start from 
divergence free initial conditions. We employ a second-order Runge-Kutta (RK2) method to 
advance equations (33) and (34) in time. We can exploit the advantage of using RK2 in the 
following two aspects. First, an integral factor can be easily introduced with the transform 

U m (t) = v m (k, t)e~» k2t , B m (t) = S m (k, t)e' XkH . (36) 

Second, aliasing errors can be reduced by introducing positive and negative random phase 
shifts at the first and second stages of RK2, respectively (Machiels & Deville, 1998). 
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The forcing term used in our simulation is the sum of two forcing functions. In Fourier 
space, it has the form 



That is, the force works only within the shell Si : 0.5 < |k| < 1.5. Here T c is a forcing term 
that is similar to, but not exactly the same as, the one adopted by Chen et al. (1993). It 
is calculated by multiplying the velocity components within shell Si : 0.5 < |k| < 1.5 by 
a factor, 7 > 1, so that before a new step of integration starts, the kinetic energy density 
within this shell is reset to E\ = 0.24. Phases of the velocity components within the shell 
are not changed. This forcing is equivalent as lengthening the velocity vectors within shell 
Si by a factor 7 — 1. Denote the increment of a velocity vector under force jF c in Fourier 
space as 5v = R + ii, where R, I are the real and imaginary parts of 5XJ. In our simulation, 
we need to inject kinetic helicity, v • V x v, into the turbulence. To do this with P C) we 
tune the angle between R and I so that they remain perpendicular to each other. Because 
kinetic helicity at k can be calculated as H(k) = 2k ■ R x I, by doing such "lengthening" 
and "angle-twisting" , we inject not only the kinetic energy but also the kinetic helicity into 
the turbulence. 

The forcing term T c maintains the energy level at the forcing scale so that the 
fluctuation in the energy development history can be small; therefore, the growth stages of 
both the kinetic energy and the magnetic energy can be studied carefully and accurately. 
However, this force does not introduce random phases into the velocity field. To be more 
realistic about the forcing in our simulation, we also add the forcing term jF b , derived 
from the forcing function used by Brandenburg (2000), as a secondary forcing function to 
introduce random phases into the velocity field. jF b has the form 



^(0.5 < |k| < 1.5) = F c + T h . 



(37) 




- k x (k x e) - i|k|(k x e) 




(38) 
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Here JF < 1 is a factor adjusted at each time step so that the kinetic energy density 
within shell Si fluctuates within ±5% of E 1 . e is an arbitrary unit vector in Fourier space. 
4>(t) is a random phase. Note that JF b (k) = jF fe (— k) so it is real, and it is helical in that 
T\> • V x jF b = —kj^ < 0, i.e., it has maximum helicity. Because jF b is tuned in such 
way that it only contributes to ±5% of the kinetic energy at the forcing scale, jF b can be 
considered as a perturbation to T c . Therefore, the advantage of using (37) as the forcing 
term is three fold: to avoid strong fluctuations of kinetic and magnetic energy density with 
time, to introduce random phases to the velocity field, and to maintain the kinetic helicity 
at certain level. 

For each of the simulation runs listed in Table 1, the initial conditions are set in the 
following way. The initial velocity field for each run is a fully developed, pure hydrodynamic, 
helical, turbulent field, obtained by applying exactly the same forcing function as that 
in Chen et al. (1993) to the Navier-Stokes equation (i.e., no magnetic field present) with 
v of that run. It is obtained and maintained helical by applying the "angle-twisting" 
method mentioned above for jF c . The initial magnetic field for each run is set up as a 
large-scale magnetic field, B, along y-direction. B is constant in both space and time. A 
magnet ohydro dynamic turbulence simulation run is then started under the forcing (37) 
with a set of initial values of the pure hydrodynamic turbulent velocity field, B and v = A, 
some of which are shown in Table 1. 

3.2. Simulation Results 

To simplify matters, we set v = A in all our simulation runs, i.e., the magnetic Prandtl 
number is 1. v{= A) and B are taken as free parameters of the numerical simulation 
and set up as initial conditions. All our simulation runs are performed on a (64) 3 spatial 
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resolution. As the MHD turbulence reaches steady state, we calculate the Reynolds 
number of the turbulence with the formula R e (= R m ) = (L)v rms /u, where (L) is the 
integral length scale of the system. If E v k is the kinetic energy spectrum, L is calculated 
as (L) = Y*kk~ l E v kl 'Y^kEvk- At t — 0, we impose a B = By into a fully developed 
hydrodynamic turbulence, and follow the MHD turbulence thereafter. We calculate the a 
coefficient through a = (v x h) y /B. Table 1 lists all the simulation runs that we obtained. 
In Figure 1, we plot the evolution of kinetic energy and magnetic energy. During the 
kinematic phase of the development, for the first few eddy turnover times, magnetic energy 
density grows as ~ t 2 , followed by an exponential growth. The growing magnetic energy 
density will impose Lorentz force on the velocity field. The kinematic phase ends when the 
Lorentz force is strong enough to significantly change the velocity field, and magnetic energy 
growth slows down. There is then a dynamic phase during which the magnetic energy 
and the kinetic energy oscillate around a certain level. During this phase, we can estimate 
the rms values of both the velocity field and the magnetic field. They are calculated 
using the formulae v rms = J2E v /3 and b rms = y2E^/3, where E v and E b are the kinetic 
energy density and magnetic energy density. In Figure 2, we plot the kinetic and magnetic 
energy spectra for the case of B = 0.1. At scales of k < 5, the kinetic energy density 
surpasses the magnetic energy density, showing that near the outer scales, the turbulence 
is largely hydrodynamic in nature. For k > 5, the kinetic energy density is smaller than 
the magnetic energy density by a factor less than threeP]. We also plot the spectra of the 

2 Exact equipartition between the kinetic energy and the magnetic energy was not found 
in our simulation. This may be due to the fact that our numerical code has a relatively low 
spatial resolution (64 3 ) so that the MHD turbulence inertial range resolved in our simulation 
is not very long. Note that the theoretical implication of the equipartition (see Blackman & 
Field, 2000) between the kinetic energy and the magnetic energy only applies to the inertial 
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absolute values of the kinetic helicity spectrum K k , the current helicity spectrum C k and 
the magnetic helicity spectrum M k (see the caption of Figure 2 for detailed definition of 
these spectra). In our simulation, we always force maximally (negative) helical flow within 
shell Si : 0.5 < |k| < 1.5. This explains the relation 2E vk (k = 1) = \K k (k = 1)| in the 
plot. For k > 1, the flow is not maximally helical. Rather, K k decreases as k increases in a 
way similar to E vk , as it decays into small scales (large k). The current helicity, C k , is also 
concentrated at large scales, k ~ 1, and decreases as k increases. In all of our simulation 
runs, we find that 90% of each of the kinetic helicity, the current helicity and the magnetic 
helicity are concentrated near the outer scales of the turbulence, i.e., k < 4. In deriving 
relations (23) and (27), we assumed that magnetic helicity and current helicity are both 
concentrated near the outer scales of the turbulence, and approximated the total magnetic 
helicity and the total current helicity with 2Lb 2 rms and \? rma jL, respectively. From our 
simulation we find that such assumptions are justified, and adopt L ~ (L). Our assumption 
of v1 ms ~ b 2 rma is not valid for all of the cases. Therefore, when comparing our estimation of 
the suppressed dynamo a coefficient (equation (27)) with the result of numerical works (for 
example, Cattaneo & Hughes (1996)), we must take this factor into account. 

To study the relation between the dynamo a-effect and the dynamics of magnetic 
helicity, we re- write equation (17) in the form 

2(v x b) -B + 2A(b- V x b) + |-(a-b) = 0. (39) 

We calculated the numerical results of each of the three quantities and the sum of them, and 
plotted them in Figure 3. Panel (a) of Figure 3 is the temporal evolution of the quantity 

range, not the dissipation range, whereas the non-equipartition found in Figure 2 and other 
works (see Brandenburg 2000, Fig. 11, and Cho & Vishniac, 2000, Fig. 7) appears mostly 
in the dissipation range. 
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2(v x b) • B. The evolution of this quantity can be separated into two stages. For t < 11, it 
first increases from to a peak value of 0.28, then decreases until it changes sign at t ~ 11. 
After t ~ 11, it oscillates around an averaged value of 0.0027 ± 0.0047. The second term 
in equation (39), 2A(b • V x b), is plotted in panel (b) of Figure 3 with A = 0.0133. This 
term represents the dissipation effect. It is a less dominant effect than the dynamo a-effect, 
which is represented by the first term in equation (39) and plotted in panel (a). The third 
term in equation (39), <9 t (a • b), is plotted in panel (c) of Figure 3. Its temporal behavior is 
similar to 2(v x b) • B. The sum of all these three quantities should be zero for our closed 
system with periodic boundary condition, and this is shown in the bottom panel of Figure 
3. 

We calculated the a coefficient with a(t) = (v x b) ■ B/B 2 and plotted it in the 
top panel of Figure 4. The evolution of the small-scale magnetic helicity density, (a • b), 
is shown in the middle panel of Figure 4. The initial value of the small-scale magnetic 
helicity is assumed to be zero. After the start of simulation, negative small-scale magnetic 
helicity is built up by the ct-effect with positive a-coefficient. The speed of this building-up 
process achieves its maximal value when the a-coefficient reaches its peak value. After 
that, the build up of negative small-scale magnetic helicity slows down as the a-coefficient 
decreases. As the positive a-coefficient approaches zero, the second term in equation (39), 
will dominate the a effect term to affect the dynamics of magnetic helicity. Panel (b) of 
Figure 4 and panel (b) of Figure 3 together show that when the dissipation term becomes 
important, the negative magnetic helicity decays. Our estimation of a(t < t c ) is 0.12 ±0.06, 
and that of a(t > t c ) is 0.008 ± 0.011. This clearly shows that the dynamo a-effect is a 
time-dependent quantity, and the constraint of magnetic helicity does not take effect on a 
until the building-up process of magnetic helicity is almost finished. 

To test how close the estimates by the models of a±, «2 and a% to the measured 
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a-coefficient from our numerical simulation, we take R = R m (^B /v^ ms J as a variable, and 

plot the measured a and the estimates of a from different models against this quantity. 
The results are given in Figure 5. In panel (a) of Figure 5, we plot a(t < t c ), a(t > t c ), oti, 
a 2 and a 3 vs. R. It is clear that a(t < t c ) is close to oti for all the values of R considered 
in our simulation runs, while 0:2 and 0:3 underestimated a(t < t c ) for large values of R. For 
R < 1, «i and a 2 give similar results. Panel (a) of this figure also shows that a 2 and a 3 
give much better estimates of the a-coefficient for t > t c , when the constraint of magnetic 
helicity on the dynamo a-effect finally enters. Panel (b) of Figure 5 shows a clear linear 
correlation between the measured critical time t c and our estimate of this quantity using 
teddy {prms/B\ , which we discussed previously in the Introduction section. 



The appearance of the electromotive force term in the equations for small-scale and 

large-scale magnetic helicity, i.e., (8) and (11), shows that the dynamo a-effect is related to 

the dynamics of magnetic helicity. But because the dynamo effect is not only determined by 

the induction equation but also the momentum equation of the velocity field, the external 

forcing term in the momentum equation will provide extra degrees of freedom to the dynamo 

a-effect, so that it is not completely determined by the dynamics of magnetic helicity. 

Instead, as we have shown in previous sections, the dynamo a-effect is largely controlled 

by the velocity field during the stage that small-scale magnetic helicity of appropriate sign 

is being pumped into large scales. At the same time as the small-scale magnetic helicity 

is being pumped to large scales at this stage, the small-scale magnetic helicity of opposite 

sign will be built up. For closed astrophysical systems of very small magnetic diffusivity, 

2 

such pumping process is in fact controlled by the dynamo a-effect, multiplied by B . If the 
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initial large scale field B is very weak, the critical time t c = t ec id y \ b rms l 'B J is initially 
very long, so the dynamical value of the a-effect, i.e., a\ of PFL and FBC, applies, and 
considerable amplification of B takes place. 

The Galaxy provides an interesting example. There, t e dd y ~ 10 7 years, b rms ~ 3 x 10~ 6 
Gauss, and B has been estimated as 1CT 13 Gauss (Field 1994). Hence initially t c ~ 10 22 
years, and dynamo action is not significantly affected by the build up of magnetic helicity. 

During this period, B will increase exponentially due to the interaction of the a effect 
and the u effect associated with the differential rotation of the Galaxy. To avoid the 
complications of the cu effect here, we will consider the simpler case of an a 2 dynamo, whose 
e-folding time is r = \/2ira. With a wavelength of lKpc and an a = 3 x 10 4 cm/sec (Field 
1994), r = 5 x 10 s years, or 50 t eddy . Presumably growth will continue until r matches the 
ever decreasing value of t c . This will occur when B = b rms /50 2 = 5 x 1CT 7 Gauss, after 
which further growth will be inhibited by the helicity constraint. To reach this stage will 
take 15 growth times, or 7.5 x 10 9 years. Thus, during a relatively long period, dynamo 
growth can occur unconstrained by magnetic helicity, and during this period, can approach 
equipartition within an order of magnitude. 

Bear in mind that this example is oversimplified. In particular, it does not take into 
account that helicity may escape through the boundaries of the system (Blackman & Field 
2000). In this case, the boundary terms f and F in equations (14) and (15) must be taken 
into account, and so one must be cautious about applying conclusions from simple models 
like we have discussed to real astrophysical systems. 

The constraint of magnetic helicity on the dynamo a-effect, as discussed by Gruzinov 
& Diamond (1994, 1995, 1996), Cattaneo & Hughes (1996) and Seehafer (1994, 1995), 
plays important role in controlling the amount of magnetic helicity pumped by the dynamo 
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a-effect from small scales to large scales. It is found in this paper that such constraining 
effect takes place after the magnetic helicity development reaches a critical time, t c . Before 
t c , we believe that the source of magnetic helicity being pumped from small scales to large 
scales may not be restricted to the Ohmic dissipation at small scales. Rather, it may come 
from the twisting, folding and stretching of the magnetic field lines by the velocity field at 
different scales. Such interactions between v and b alter the topology of the magnetic field 
in such a way that net magnetic helicity, which represents the magnetic field line topology 
(Moffatt & Tsinober 1992), can be built up. Another source of magnetic helicity can be 
from outside of the system if there are open boundaries. During such non- stationary stage, 
induction equation (6) alone cannot give a complete picture of the dynamo a-effect, and 
the dynamical studies by PFL and FBC may provide a valid a-coefncient, cci, which cannot 
be equated to a 2 - When the building up of magnetic helicity approached its upper limit set 
by the realizability condition (22), the dynamical pumping effect described by a.\ starts to 
be constrained by the magnetic helicity conservation. Eventually, if the system is closed, 
cti alone (in the form of relation (2)) is not enough for a complete picture of the dynamo 
a-effect, and a 2 must be introduced to get a^. 

Another motivation of our work is the numerical study of magnetic helicity by Stribling 
et al. (1994). In their simulation, they find that the electromotive force is not suppressed 
for the first few eddy turnover times in a 3D MHD turbulence with an imposed moderately 
strong large-scale magnetic field and an R m — > oo. Our work is an extension of the work 
by Stribling et al. in the following aspects: first, we introduced a critical time to separate 
the non-suppressed stage from the suppressed stage of the a-effect; second, we numerically 
studied the dependence of the a-coefficient on the values of B, b rms , r e ddy and other 
quantities of the MHD turbulence. By doing so, we argued that the time behavior of the 
magnetic helicity dynamics plays an important role in the dynamo a-effect; therefore, one 
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cannot simply ignore the <9t-terms when applying the magnetic helicity equations, (17) and 
(18), to real astrophysical systems. 



We studied the constraint of magnetic helicity on the dynamo a-effect with 3D direct 
numerical simulation under periodic boundary conditions. The dynamics of magnetic 
helicity affects the dynamo a-effect only after the magnetic helicity at small scales is built 
up and the magnetic helicity dynamics enters a stationary state. Such building-up process 
can be understood as a pumping effect of the dynamo a-effect, and the a-coefficient during 
this non-stationary pumping stage can be estimated according to the model by Pouquet et 
al. (1976) or Field et al. (1999). As the small-scale magnetic helicity is built up to the level 
limited by the realizability condition, the a-effect is quenched, as suggested by Gruzinov & 
Diamond (1994, 1995, 1996) and Cattaneo & Hughes (1996). The a-coefficient during such 
magnetic helicity constraining stage can be estimated according to the model by Seehafer 
(1994, 1995, see also Blackman & Field 2000). A critical time, t c ~ t eddy (b rms /B\ , is 



introduced to separate these two stages. 

We benefited from our discussions with E.G. Blackman, A. Brandenburg, B. Chandran 
and R. Kulsrud. We thank an anonymous referee for insightful comments. 
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Table 1. Measurements of various physical quantities for different simulation runs 





Run I 


Run 11 


1 > TTT 

Kun 111 


Run IV 


Run V 


B 


0.316 


0.1 


0.0316 


0.0316 


0.224 


«/(= A) 


0.0133 


0.0133 


0.0071 


0.0133 


0.0133 


v rrns 


0.473 ±0.006 


0.498 ±0.008 


0.50 ±0.01 


0.487 ±0.005 


0.489 ±0.002 


forms 


0.35 ±0.04 


0.26 ±0.04 


0.24 ±0.07 


0.19 ±0.05 


0.27 ±0.07 


Rrn 


29.8 ±0.9 


31.2 ± 1.4 


54.9 ±2.1 


28.5 ±0.6 


29.6 ±0.9 


teddy 


1.6 ±0.1 


1.6 ±0.2 


1.6 ±0.7 


1.7 ±0.4 


1.7 ±0.2 


Tcor 


0.64 ±0.08 


0.8 ±0.3 


0.7 ±0.4 


0.5 ±0.2 


1.0 ±0.4 


t c 


~ 11 


- 32 


- 85 


- 120 


- 20 


a(t < t c y 


0.12 ±0.06 


0.21 ±0.08 


0.26 ±0.09 


0.24 ±0.26 


0.13 ±0.09 


a(t > t c y 


0.008 ±0.011 


0.10 ±0.11 


0.20 ±0.09 


0.17 ±0.22 


0.04 ±0.02 




0.11 ±0.05 


0.23 ±0.05 


0.16 ±0.12 


0.22 ±0.05 


0.16 ±0.0.01 


a 2 h 


0.006 ±0.006 


0.11 ±0.07 


0.14 ±0.05 


0.18 ±0.17 


0.01 ±0.01 


a 3 h 


0.007 ±0.003 


0.09 ±0.03 


0.13 ±0.04 


0.19 ±0.06 


0.022 ±0.002 



a Tcmporal average of the numerical a-coefficient (= (v x h) y /B) before and after the 
critical time t c . 

b Estimates of the a-coefficient by three different models (see equations (37), (38) and 
(39) for details). 
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Fig. I. — Temporal evolution of the kinetic energy density (thin solid line) and magnetic 
energy density(thick solid line) for a few run cases. 
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jj dk'|C k (k')|/P 2 dk'|C k (k')| 
jj dk'|M k (k')|/lf dk'|M k (k') 
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Fig. 2. — Top panel: kinetic energy spectrum, E vk (k) = \ YX-k\ l^(^')| 2 ; magnetic energy 
spectrum, E bk (k) = ^Y^k-o.l IK^OI 2 ! absolute value of the kinetic helicity spectrum, where 
K k (k) = X^oJik' ' (v(k') x v*(k')); absolute value of current helicity spectrum, where 
Ck(k) = X^toJik' • (b(k') x b*(k')); absolute value of magnetic helicity spectrum, where 
M k (k) = Ylk-o 5(^(^0 • b*(k')) and a is the vector potential of b; kinetic helicity spectrum 
of maximally helical flow, 2kE vk (k). Relative error for E vk {k) ranges from ±5% to ±40% for 
different fc's, with k — 1 has the smallest error. Relative error for E bk {k) ranges from ±30% to 
±50%. Bottom panel: cumulative spectra of E vk) E bk , \K k \,\C k \ and \M k \. Data are collected 
from Run II, averaged from time t = 52.5 to t — 100.5. For this case, (L) = 0.80 ± 0.02, 
v rms = 0.498 ± 0.008; therefore Reynolds number is R e = R m = 29.8 ± 0.9. 
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Fig. 3. — Temporal evolution of the following quantities: (a) 2(v x b) • B; (b) 2A(b • V x b); 
(c) d t (a ■ b); (d) sum of above three quantities. Data are collected from Run I. 
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Fig. 4. — Panel (a): ct(t) = ^ V ^' B . Panel (b): magnetic helicity density (a - b). Panel (c): 
— Jq 2a (r)B 2 dr (thick solid line) and —2 Joi](h ■ V x b)(thin solid line). Panel (d): negative 
kinetic helicity density, — (v- V x v) (thick solid line); current helicity density, (b-Vx b)(thin 
solid line). t c ~ 11. Data are collected from Run I. 
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Fig. 5. — Panel (a): The measured a(t < t c ), a(t > t c ), and estimates of a based on different 
theoretical models (see text for details). Panel (b): the critical time t c measured from the 
simulation vs. the estimate based on our model (see section 1 of main text). The correction 
factors to t c are estimated by reading the spectrum of magnetic helicity and calculating 
the ratio v rms / 'olfbc- The correction factors are C = 1.5,2, 1,2.5 and 3 for Runs I to V. 
Error bars shown in panel (b) are the standard deviation of our calculations for t c . For 
the statistical standard deviation calculations of other quantities in this figure, the reader is 
referred to Table I for more details. 



